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16. Abstract 

The importance of a uniform current density profile i" ®^5pm^lifetirnf ^A ?IlidS?e^t1me^°" 
thruster is discussed in terms of thrust level and accelerator system lifetime. A residence 
aooroach is used to explain the nonuniform beam current density profile of the SERT 

fipiH thruster and to propose a magnetic field modification (the cusped magnetic field 

having a highly uniform beam profile, good performance, and a tu„ thruster 

Exoprimental results are presented for a wide range of the discharge chamber length. The thruster 

desianed for this investigation was operated with a cusped magnetic field as well as a 
aeomLrv and the cusped field geometry is shown to be superior from the standpoint of beam profile 
uniformity performance, and double ion population. Beam profile measurements 6 mm downstream of the 
accelerator grid indicate the beam profile flatness parameter (which 

uniformity of the profile) is 0.7, and throttled thruster operation results in a beam flatness 
parameter of 0 82.^ At about the same thrust and performance levels, the beam flatness parameter of 
thJ cSSed fie d thruster is 40 percent higher than the SERT II thruster value. This implies a 
SS percent increase i^the accelerator system lifetime due to a reduction in localized erosion of the 
accelerator grid caused by the impingement of charge-exchange ions. 
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I. INTRODUCTION 


Propelling space vehicles by electrical means was proposed as early 
as 1906 by Robert H. Goddard [1], and since that time numerous electric 
space propulsion concepts have been introduced [1,2]. Among these pro- 
posals the electrostatic acceleration of positive ions is generally 
considered the most promising for high-impulse space missions. In 
particular, the electron-bombardment ion thruster introduced in 1960 
[3,4] remains today the most likely choice for primary and many auxil- 
iary space propulsion needs [5]. A discussion of the historical 
development and current technology status of electron-bombardment ion 
thrusters can be found in the literature [5,8]. 

Thruster Lifetime 

The low- thrust nature of electron-bombardment ion thrusters 
requires long mission times to achieve the desired spacecraft velocity 
increment. This requirement dictates thruster designs which can operate 
for extended periods of time; perhaps as long as 20,000 hours [9,10]. 
Designing thrusters to operate for such long mission times poses many 
unique and challenging problems for the designer and requires life- 
testing of critical thruster components under simulated space conditions. 
The operational lifetime of a thruster is usually determined by the 
failure of one of its components rather than failure of the thruster as 
a whole. However, in some instances an early failure may be due to an 


in l,or<K. t ion bol.wiHni Iwo or mote comiiononl:'- v^Iik.I' mu', i v : ei i m / Ii.kI iiiikJi 
lon<jt.T 1 i tciLiiiio',) . lor o;<.iiiiiilis I'olii Ihrns m. >. ■ .1 
soarnt.rii ft , which was lannoliC'''.l in 19/lJ 111 It IciiOi' •i‘> <t tc.uli "i 
olectricai shorts in the accelerator system. Ihe i aose of Liu* failuM; 
was due not to poor accelerator system desiyn bu i. to toe close pruxiim^ 
of the accelerator yrid and the neutralizer catiiodc’ assembly. Lovr- 
eneryy ions created by the neutralizer discliaryo vtere accelerated to tm 
downstream side of the accelerator grid resulting in localized sputter- 
ing. Metallic fragments produced by the erosion process eventually 
became lodged between the screen and accelerator or ids resulting in an 
electrical short. The problem has since been eliminated simply by 
relocating the neutralizer cathode. 

The extension of the electron-bombardment thruster lifetime from 
the 100 or so hours of the earliest designs to around 10,000 hours for 
the more recent designs [12j has involved a somewhat systematic process 
of prolonging the life of individual components. Diiring a development 
program the same thruster component may emerge periodically as the 
lifetime limitation. An example is the accelerator system which, as 
discussed earlier, failed and caused the termination of the SERT 11 
mission. The source of this failure was easily conjectured from the 
results of subseguent ground tests, and a simple solution was identified 
and incorporated into current thruster designs. Tliere are, however, 
additional accelerator system erosion phenomena which arc recognized 
today as posing potential life limitations: (I) erosion of the accel- 

erator (jrid caused by chanie extiiange impiiKiement , {?.) accel er<M or grid 
('rosion caused by direc t ion impin'jement , and (.1) erosion of the sc r(?en 
fi r('Mi I i ')f l)nm|hir(liii('nl l*v dnuM y ( inri',. 
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Charge-exchange ions are produced when an ion and neutral passing 

through the accelerator system exchange identity. This process occurs 

when an electron leaves the slow neutral and attaches itself to the 

★ 

fast ion. The mercury charge-exchange cross section is about 
-19 2 

6 X 10 m for the ion energy range of interest to electric propulsion, 
and this is by far the largest of any of the cross sections associated 
with the reactions occurring during thruster operation. The products 
of the charge-exchange reaction are a fast neutral and slow ion, and if 
the reaction occurs near the negative accelerator grid the low-energy 
ion will be accelerated toward the grid and upon impact may cause con- 
siderable sputtering damage. The severity of the problem is illustrated 
in Figure 1 which shows the charge-exchange erosion pattern observed on 
the downstream side of an accelerator grid after having been operated 
for about 3,000 hours. Deep triangular pits are formed in the regions 
surrounded by three adjacent holes where the impinging current density 
is highest. Although these deep pits are first to erode completely 
through, failure of the accelerator grid results from the eventual 
breakthrough of the shallower troughs in the regions adjacent to two 
holes. The upstream side of the accelerator grid also undergoes charge- 
exchange erosion; but at a much lower rate. 

Direct ion impingement upon the accelerator grid (as opposed to 
charge-exchange ion impingement) may result during thruster operation at 
high plasma density conditions. Under these conditions the plasma 
sheath tends to penetrate the screen grid holes and assumes a curvature 
★ 

Mercury is generally chosen as a propellant because of its low 
ionization potential, high molecular weight, and ease in handling and 
storing. 








5 


wh’ch can direct plasma ions into trajectories torminatinq on the 
accelerator grid. Low plasma density operation can also cause direct 
ion ifiipingement as a result of crossover trajectories. Crossover occurs 
when the plasma sheath recedes away from the screen grid aperture and 
assumes a highly concave geometry which tends to focus ions along tra- 
jectories which cross over the axis of the hole. The low-density region 
near the chamber wall can also cause direct ion impingement when the 
plasma sheath moves upstream to compensate for the reduction in plasma 
density and assumes a non symmetrical shape which directs ions onto one 
side of the accelerator grid hole. Due to the nonuniform plasma 
density profile which exists inside most discharge chambers, the first 
effect tends to occur in the center of the chamber where the plasma 
density is high while the other two effects occur near the edge of the 
chamber where the plasma density is lower. 

Screen grid erosion due to the impingement of doubly charged ions 
may also represent a lifetime limit on accelerator systems. Although 
the ion energy is considerably less than in the case of charge exchange 
and direct impingement, the amount of material available for erosion is 
substantially reduced also. This is because good thruster performance 
requires thin screen grids to minimize recombination losses and thick 
accelerator grids to minimize neutral particle loss. The long-term 
effects of double ion erosion can be serious [13], and perhaps the best 
approach to eliminating this problem is to reduce the number, as well 
as the energy, of the double ions present in the discharge chamber. 
Substantial progress in this area has been realized recently by use of 
small-hole accelerator grids [14]. These electrodes have accelerator 


holes which are only about ono-half the diameter of the screen grid 
hole and, at the same propellant flow rate, result in a higher neutral 
density within the discharge chamber. This allows efficient operation 
at much lower discharge voltage levels which results in a reduction in 
double ion production rate and particle energy. 

Present Investigation 

The specific objective of the present investigation is to increase 
the accelerator grid lifetime by reducing local charge-exchange ion 
erosion. The most direct approach to the problem is to reduce the rate 
of production of charge-exchange ions. Some insight into how this 
might be accomplished can be gained by considering the expression for 
the accelerator grid charge-exchange current 

^CE "" ^o'^b'^Ce'^ 

where = accelerator grid charge-exchange current 
n^ = neutral particle density 
jg = beam current density 
Ogg = charge-exchange cross section 

V = charge-exchange ion production volume 

The charge-exchange cross section is determined by the choice 
of propellant and the optimum specific impulse for the mission, while 
the volume V is determined by the accelerator system design parameters 
(grid spacing and thickness, location of the neutralization plane, and 
beamlot diameter). Therefore, for a given accelerator system the only 
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other variables remaining which are at the designers disposal are the 
neutral density and beam current density. For a given propellant flow 
rate the beam current density is proportional to the propellant utiliza 
tion* Hy. and the neutral density is proportional to the quantity 
1 - riy* The charge-exchange current can thus be expressed as a simple 
proportional ity 


'CE ” %<' - 

From this expression one can see that the charge-exchange current is 
zero when the propellant utilization is either zero or unity. The 
former case corresponds to all the propellant exiting the thruster in 
the unionized state and is of no interest to space propulsion. The 
latter case is unrealistic from physical considerations since a pro- 
hibitively large power supply would be required to ionize 100 percent 
of the propellant. This can be illustrated by considering a typical 
thruster performance curve presented in Figure 2 . The performance 
parameter plotted as the ordinate is the discharge power loss, or the 
energy expended in producing a beam ion, and is calculated by use of 
the expression 


P 


I V 
anode anode 


beam 


(3) 


^Propellant utilization is defined as the fraction of the total 
propellant flow rate which exits the thruster in the ionized state. 
This parameter is of great importance in the electric propulsion field 
since particles leaving the accelerator system as neutrals p-^oduce 
negligible thrust. 





where P = discharge power loss, eV/ion 

'anode ' 

"anode ” potential 

'beam ' 

The bend in the performance curve is referred to as the "knee" of the 
curve and generally represents the optimum operating point. Operation 
to the right of the performance curve knee requires a substantial 
increase in power supply weight with a relatively small reduction in 
propellant weight. Operation to the left of the knee results in a 
substantial increase in propellant mass accompanied by only a slight 
decrease in power supp''y mass. Equation 2 indicates a higher charge- 
exchange current will generally result from thruster operation to the 
left of the knee. 

At first glance it appears little can be done to reduce the 
charge-exchange ion production. However, closer examination of Figure 
1 indicates the erosion damage is strongly dependent on the radial 
coordinate, and the erosion rate is greatest on the thruster centerline. 
Current density measurements in the thruster exhaust beam reveal the 
profile is quite peaked, and the maximum current density occurs on the 
beam centerline. These observations, along with the dependence of the 
charge-exchange current on the local beam current density as illus- 
trated by Equation 1, explain the radial variation of the charge- 
exchange erosion. A quantity which describes the nonuniformity of the 
beam current density profile is the flatness parameter F which is 
defined as the ratio of average-to-maximum beam current density. 
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Assuming the beam is axisymmetric the flatness parameter can be 
calculated by the relation 

R 

2irrj(r)dr 


where F = beam flatness parameter 
r = radial coordinate 
R = beam radius 
j = current density 

Figure 3 presents an example of two beam piofiles which have the 
same average current density. The uniform profile has a flatness 
parameter of 1.0, and the nonuniform profile (a cosine distribution) 
has a value of F equal to 0.46. This nonuniform profile is typical of 
state-of-the-art divergent field chamber designs and indicates the 
current density at the beam center is more than twice the value that 
would exist with a uniform profile having the same integrated current. 

The dependence of the accelerator grid erosion rate on local beam 
current density suggests the grid (assumed to have homogeneous compo- 
sition and thickness) will wear through first in regions where the 
current density is a maximum. Experimental measurements indicate the 
current density is a maximum on the beam centerline, and, thus, the 
wear rate pattern illustrated in Figure 1 is to be expected. These 
observations suggest an improvement in accelerator grid life could be 
realized by distributing the charge-exchange ion erosion more evenly 





"T 
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across the grid. The correlation between erosion and current density 
profiles suggests a more uniform erosion rate could be realized by 
increasing the beam flatness parameter. 

An additional incentive for achieving a uniform current density 
profile can be illustrated by use of the thrust equation for an 
electrostatic accelerator 


T = j._A(2mvVq)^^^ 


'ave 


( 5 ) 


where T = thrust 

i = average current density 
'^ave 

A = beam area 

= net accelerating potential 
q/m = propellant charge-to-mass ratio 

For a given propellant, exhaust velocity, and thruster diameter, the 
thrust equation reduces to a simple proportionality between thrust and 
average beam current density. By introducing the beam flatness 
parameter F, the proportionality can be written as 


Fj 


max 


( 6 ) 


Since the accelerator grid lifetime is determined by the maximum beam 
current density. Equation 6 indicates the thrust is directly propor- 
tional to beam flatness for a given accelerator grid lifetime. Uniform 
beam profiles are therefore seen to be highly desirable from the 
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standpoint of (1) increasing accelerator grid lifetime at a given 
thrust level by reducing localized charge-exchange ion erosion or (2) 
increasing the thrust level for a given accelerator grid lifetime. 
Either improvement is of interest to the thruster designer, and 
together they form the basis of this investigation. 

Unless otherwise noted, SI (rationalized mks) units are used 
throughout this paper. 







application of a neutral residence time theory to a modern discharge 
chamber design and suggests a magnetic field modification which should 
produce a more uniform ion beam current density profile. The final 
sections discuss the effects this modification is expected to produce 
on thruster performance and the discharge chamber double ion content. 


Thruster Operation 

(Contemporary electron-bombardment ion thrusters consist of a 
cylindrical discharge chamber with a set of grids mounted at the down- 
stream end as illustrated in Figure 4. The cylindrical chamber wall 
serves as the anode, and the upstream end of the chamber houses a 
hollow cathode [7]. The cathode is mounted within a soft-iron pole 
piece and is shielded from the bulk of the thruster discharge by use of 
a circular baffle. This divides the thruster into two discharge regions 
and provides sputter erosion protection for the hollow cathode. 


A 


A 


MAGNET 


PROPELLANT 



Figure 4. Electron-bombardment ion thruster. 



16 




Operational lifetimes on the order of 10,000 hours have been realised 
with this arrangement [15]. 

The ion accelerating system consists of two closely-spaced grids 
containing many small diameter holes. The grids have a slightly 
spherical shape which gives the desired stability under the thermal 
loading produced by temperature gradients in both the radial and axial 
directions. The screen grid has as large an open-area fraction as 
structural and fabrication considerations will permit. Usually this 
amounts to about a 70 percent open-area fraction which is sufficient to 
extract nearly all the ions reaching the accelerator system. In order 
to minimize ion recombination on the walls of the cylindrical holes, the 
screen grid is made as thin as structural and fabrication considerations 
will allow. The accelerator grid, on the other hand, usually has a low 
open-area fraction and can be much thicker, both of which help to 
minimize the loss of neutral particles. The extracted ion beam is 
current and space charge neutralized to eliminate spacecraft charging 
and space charge accumulation within the ion beam. Neutralization is 
achieved by the use of a hollow cathode neutralizer which is similar 
in design to the hollow cathode used inside the thruster. 

During typical thruster operation with mercury propellant, up to 
90-95 percent of the neutral gas is injected into the thruster at the 
upstream end of the discharge chamber. The remainder of the propellant 
enters through the hollow cathode. The neutral gas density near the 
accelerator system is of the order of 10 ' and is essentially 
independent of propellant flow rate at efficient operating conditions 
[6]. The low neutral density results in free-molecular flow of the 
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a I in, c species, and the neutrals acquire the mean chamber wall 
temperature of around 500 °K. 

Ionization of the neutral propellant atoms, which reside in the 
ground energy state and various excited states, is accomplished by 
electron impact involving two groups of ionizing electrons. The first 
group is comprised of the electrons which are accelerated from the 
hollow cathode into the main discharge region, and these electrons 
acquire an energy roughly equal to the difference in the plasma poten- 
tials of the two regions. The difference in potential is usually of 
the order of 30 V, and the high-energy electrons in this group are 
referred to as monoenergetic, or primary, electrons. The second group 
of electrons originates from the inelastic interactions of the primary 
electrons with neutrals and ions. These interactions reduce the energy 
of the primary electrons and release low-energy secondary electrons. 

The low-energy electrons rapidly randomize and approximate a Maxwellian 
distribution of energies characterized by a temperature of the order of 
5 eV. The Maxwellian group usually contains 90-95 percent of the plasma 
electrons. However, this value can change with thruster size and 
geometry and is quite sensitive to thruster operating conditions. The 
remainder of the electrons are primaries which, although few in number, 
contribute a large percentage of the ionizations due to their higher 
energy. 

During operation the propellant gas is only partially ionized, and 

17 .7 

an electron density of about 10 m is typical. The relatively low 

I electron temperature and high density results in a Debye length of 

_ti 4 

5 X 10 m. The Debye number is about 6 x 10 , and under these 

conditions the plasma may be considered quasi-neutral. 
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Magnetic Field 

The mean free path for electron-atom ionizing collisions is about 

1 8 " 3 

20 m for 30-40 eV electrons and a neutral density of 10 m . Since 
the path length is much greater than the thruster dimensions, a magnetic 
field is employed to prevent electrons from escaping to the anode with- 
out first encountering a neutral. Magnetic field strengths of about 
5 X 10'^ W/m^ are typical of most electron-bombardment thrusters and 

_3 

result in an electron gyro radius of about 4 x 10 m for 30 eV 
electrons. The ions, on the other hand, are essentially unaffected by 
the magnetic field due to their large mass. 

The shape of the ion beam current density profile is highly in- 
fluenced by the magnetic field geometry which exists inside the dis- 
charge chamber. This effect is illustrated in Figure 5 which indicates 
three different discharge chamber designs and their associated beam 
current density profiles. Studies using a divergent magnetic field 
geometry [16] have indicated a good statistical correlation exists 
between the shape of the beam current density profile at the accelerator 
grid and the ion number density profile near the screen grid. Similarly, 
a good qualitative correlation between the ion density and beam current 
density profiles was observed with the radial field thruster design [17]. 

The relationship between the magnetic field geometry and ion 
density profile can be explained on the basis of electron transport 
along, and transverse to, the magnetic field lines. In the axial 
magnetic field thruster illustratec in Figure 5a, the high-energy 
electrons are emitted by a refractory-metal cathode located on the 
thruster centerline. The emitted electrons have direct access to the 
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central core of the chamber and reach the remainder of the chamber cror.<i 
section by collisions with neutrals and, to a lessor extent, ions. The 
inelastic collisions with neutrals and ions results in primary electrons 
losing much of their energy. As a result the high-energy electrons are 
concentrated in the central core region, and the ion production rate and 
ion density decrease radially towards the anode. 

The divergent magnetic field geometry illustrated in Figure 5b 
allows the energetic electrons access to the entire chamber cross 
section by traveling along the magnetic field lines which originate near 
the baffle aperture region. The primaries are no longer confined to the 
central portion of the chamber but have access to the entire bell -shaped 
region as illustrated in the figure. This region is known as the 
primary electron region, and the upstream boundary is called the 
"critical field line." The critical field line is defined as the inner- 
most field line which intercepts the anode. The basis for this defini- 
tion is simply that once an energetic electron crosses this line it is 
collected by the anode. This occurs because the potential drop in the 
anode plasma sheath is insufficient to reflect the high-energy electrons. 
As a consequence the high-energy, or primary, electrons are confined to 
the "primary electron region." The Maxwellian electron temperature also 
drops beyond the critical field line since the anode represents a sink 
for the electrons in the tail of the distribution. The low-energy 
electrons which exist outside the primary electron region are inefficient 
in ionizing neutrals, and consequently the plasma density in this region 
is much lower. Experimental results have shown that thruster performance 
is determined by the plasma conditions existing within the primary 
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eloctron region and dir.ch.irge chainlter (|(‘iiiiK l.ry (;h.m<|Os extorridl to 

ttiis region have essentially no effect on l.hrust.er |)erformanc;e (6,19j. 

A comparison of the Ixiaiii current density i^rofiles of the axial and 
divergent field thrusters indicates a substantial hiiprovoinent in the 
uniformity of the profile has been realized as a result of expanding the 
primary electron region to cover the entire chamber cross section at the 
downstream end of the thruster. 

Another discharge chamber geometry which results in a further flat- 
tening of the ion beam profile is the radial field thruster illustrated 
in Figure 5c. This geometry has an extreiDely divergent, or radial, 
magnetic field arrangement. High-energy electrons emitted by the 
cathode are free to gyrate about field lines which are directed radially, 
and the electrons diffuse axially to the anode which is located in the 
upstream end of the chamber. Plasma property measurements in this 
thruster reveal a uniform density profile in the radial direction, and 
this uniformity is reflected in the beam profile shape illustrated in 
Figure 5c. The beam flatness parameter of this profile is about 0.51 
and was the highest value reported in the literature at the time this 
investigation was initiated. 

Residence rime Theory 

A comparison of the sliape of the critical magnetic field line and 
the beam current density profile of the div(?r(|ent field thruster 
illustrated in Figure 5b reveals a striking similarity. The h(?am 
current density profile; is bel 1 -shigxi'd ; much like the upstream boundary 
of the primary electron region. This observation su(|g('-,ts the inn 
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donsity nocir the screen grid is a function of the distance between the 
critical field line and the grids; a reasonable result considering the 
probability of a neutral encountering an ionizing electron as it drifts 
through a specified region. 

If one considers neutral trajectories which are parallel to the 
thruster axis, the length of time a neutral atom resides in the primary 
electron region is given by the expression 


where r = neutral residence time 

L = distance between the critical field line and screen grid 

V = neutral particle speed 
0 

In order to encounter an ionizing collision the neutral atom must, on 
the average, reside in the primary electron region for a period equal 
to the characteristic time for ionization. The ionizing collision fre- 
quency for a neutral ground state atom interacting with primary and 
Maxwellian electrons is given by the expression 

v = naV+n<TaV^ (8 

^ p p m ^ m ^ 

where v = ionizing collision frequency 

np = primary electron density 

n^jj = Maxwellian electron density 

V = Maxwellian electron speed 
m 


Vp = primary electron speed 
fj = ionization cross section 

The symbol <Q > denotes the average of an energy-dependent quantity 
Q for the Maxwellian distribution function. Equating the neutral 
residence time with the reciprocal of the ionization frequency results 
in the following expression for the minimum depth of the primary 
electron region 


L = ? — (9 

n„jV„ + n„<aV„> 
p p m ^ m ' 

The neutral residence length L was calculated using Equation 9 and 
the measured plasma properties of a SERT II divergent field thruster. 

An iron filings map was used to determine the distance between the 
critical field line and the screen grid as a function of thruster 
radius. This distance was normalized by the calculated value of L 
for each radial position and in this form represents the probability 
for ionization. The results of these calculations are presented In 
Figure 6 along with the measured beam current density profile of this 
thruster. The similarity between these profiles suggests the neutral 
residence time consideration is a viable approach to shaping the beam 
current density profile. This implies the beam profile can be made 
more uniform by increasing the neutral residence time at the periphery 
of the discharge chamber. Since the residence time is determined by 
the distance between the critical field line and screen grid, moving 
the field line upstream at the thruster periphery should produce the 





desired effect. However, this modification will increase the volume of 
the ion production region, and the effect of this increase on thruster 
performance must also be taken into account. 

Performance Considerations 

The SERT II thruster is a high-performance engine which was 
developed during an extensive test program involving over 100 discharge 
chamber configurations [20]. The goal of the present investigation is 
to improve the beam profile flatness of this thruster as a means of 
extending the accelerator system lifetime. At the same time it. is 
desirable to maintain, or perhaps improve upon, the high performance 
level of the SERT II thruster. For this reason, expressions relating 
the thruster power loss and propellant utilization to the discharge 
chamber geometry will be developed and used to predict thruster modifi- 
cations which result in a uniform beam profile while maintaining good 
overall performance. 

An expression for the propellant utilization efficiency can be 
derived by considering the single ion production and loss rates. It is 
assumed that uniform plasma properties exist within the ion production 
region, and that single ions are produced from the atomic ground state 
as a result of primary and Maxwellian electron impact. The ions are 
lost at the boundary of the production region at a rate given by the 
modified Bohm criterion [21]. Equating the production and loss rates 
for steady-state operation results in the following expression for the 


where = ion density 

n = neutral density (ground energy state) 

0 

= modified Bohm velocity 
V = ionization collision frequency 
V/A = primary electron region volume-to-surface area ratio or 
characteristic length 

The propellant utilization expression is obtained by multiplying the 

★ 

ion flux density by the screen grid effective open area and 
dividing by the propellant flow rate m 
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n^vA.V/A 
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( 11 ) 


The neutral loss rate and propellant flow rate are related to the pro- 
pellant utilization by the expression 

Oq = (1 - ny)m C2) 


where n is the neutral loss rate and is given by the free-molecular 
0 

flow expression 


*The effective open area of the screen grid is sometimes greater 
than the geometric open area due to the plasma sheath assuming a con- 
cave shape near the upstream side of the screen grid. The effect of 
this condition is to deflect ions into the beam which would otherwise 
have recombined on the screen grid. 


Aq is the effective sharp-edged orifice area of the accelerator system. 
Combining Equations (11), (12), and (13) results in the desired expres- 
sion for the propellant utilization 
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(14) 


In deriving Equation (14), it was assumed that single ions were 

produced from the atomic ground state. However, the excited states of 

1 3 

mercury (particularly the 6 P.| resonance state and 6 P 2 metastable 
state) contribute substantially to the total production rate [22]. To 
account for this, the ground state reaction rate n^v should be replaced 
by the total reaction rate which includes the excited states. In the 
final result, the collision frequency v can be replaced by an effective 
collision frequency v which is given by the following expression 
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where the subscripts r and m refer to the resonance and metastable 


states, and the summation is over all excited states which contribute 


28 



significantly to the production of single ions. The collision fre- 
quencies and vj. are given jy Equation (8) with the cross section o 

i i 

replaced by the cross sections and o^. 

As a first approximation the discharge power, or beam ion produc- 
tion cost, can be expressed as the product of the plasma ion production 
cost C. and a factor which is the ratio of primary electron region 
surface area to screen grid effective open area [18] 

P - C.A/A. (16) 

The plasma ic. production cost represents the energy required to 
produce a plasma ion and consists primarily of inelastic collisional 
losses. 

If one assumes the plasma properties remain constant, or vary in 
such a way that C- and v remain constant, then Equations (14) and (16) 
can be used to demonstrate the effect of discharge chamber geometry on 
thruster performance. For example, Equation (16) shows the discharge 
power will remain constant if the ratio A/A. is held fixed. This re- 
quirement and Equation (14) indicate the propellant utilization also 
remains constant if the ratio A^/V is fixed. These results are useful 
scaling relationships for the thruster designer. For instance, if the 
primary electron region has a cylindrical shape the performance 
parameters given by Equations (14) and (16) can be simplified to the 
following form by introducing the area and volume relationships for a 
cylinder of length ■' and radius r 



4)^ and 4)^. are the effective open area fractions of the accelerator 
system to neutrals and ions, respectively. From these expressions one 
can see the importance of a large ion open area fraction in reducing 
discharge losses and increasing propellant utilization. The dependence 
of propellant utilization on the product Vq4)q illustrates the importance 
of low chamber wall temperature, high molecular weight, and low neutral 
open-area fraction. Both the discharge power and propellant utiliza- 
tion depend on thruster radius in such a manner that large diameter 
thrusters should perform better than small diameter ones; a well- 
established experimental result. The length, on the other hano, is 
seen to have opposing effects on the performance parameters. This 
suggests small diameter thrusters can achieve acceptable propellant 
utilization efficiencies (at the cost of higher discharge losses) by 
increasing the chamber length. This prediction is also apparent if one 
considers the variation in the length-to-diameter ratios of optimum 
thruster designs [19]. 

The results of the preceeding analysis can be used for the some- 
what different problem which is of interest to the present investigation. 
In this case the thruster radius is fixed, and the primary electron 
region can be modeled by a cylinder which lies along the axis of the 


thruster as illustrated in Figure 7a. Due to the critical magnetic 
field line geometry, the effective radius of the cylinder is consider- 
ably smaller than the chamber radius. This confinement of the high- 
energy electrons to the center of the chamber results in a nonuniform 
beam current density profile as illustrated in the figure. In order to 
achieve a more uniform beam profile, the radius of the cylindrical 
primary electron region can be increased by means of a magnetic field 
modification. The basic difference between this problem and that of 
increasing the thruster diameter is that in the present case the same 
grids are used and thus the open area for neutral loss remains constant 
as the primary electron region is expanded radially outward. The 
effect this change will have on the thruster performance parameters 
can be predicted by use of the expressions derived previously. 

Equations (14) and (16) were used to calculate thruster performance 
as a function of cylinder length for fixed values of cylinder radius. 

The quantity V^A^/Av in Equation (14) was evaluated using the SERT II 
thruster geometry and propellant utilization. Calculations were made 
for cylinder radii of 4 cm and 6 cm which were chosen on the following 
basis. The SERT II thruster primary electron region has a volurne-to- 
area ratio of 1.4 cm and the length at the thruster centerline is 9 cm. 
An equivalent cylinder having this same length and volume-to-area ratio 
has a radius of 4 cm. The 4 cm value also represents the radius at 
which the beam current density drops off rapidly as illustrated in 
Figure 6. The 6 cm value was chosen as representative of the radius to 
which the primary electron region can be expanded within a 15 cm 
diameter discharge chamber while maintaining adequate plasma confinement 
at the cylindrical chamber boundary. Different choices for the values 



of r would have no qualitative and little quantitative el feet on the 
results of the calculations. The performance parameters calculated in 
the manner described above are presented in Figure 8. The numbers 
adjacent to the symbols on these curves indicate the primary electron 
region length in cm. It is interesting to note the similarity between 
these curves and performance plots obtained in an operating thruster 
(e.g.. Figure 2). The increase in primary electron region length has 
the same effect as increasing the power input, while the increase in 
radius has the same effect as increasing the propellant flow rate. The 
data of Figure 8 indicate two important results; (1) increasing the 
diameter of the primary electron region produces substantial perform- 
ance gains and (2) an optimum length exists for each radius which 
represents a 'ompromise between discharge losses and propellant utili- 
zation. This latter result is analogous to the performance curve 
"knee" discussed earlier in the chapter. 

The results of the analysis described in this and the previous 
section suggest the modification to the SERT II discharge chamber 
geometry illustrated in Figure 7b. The changes consist of (1) a 
magnetic field modification which increases the radius of the primary 
electron region and (2) a reduction in the primary electron region 
length. The anticipated results of these modifications are an improve- 
ment in both the uniformity of the ion beam current density profile and 
the thruster performance. The improvements predicted by the simple 
model are, of course, based on the assumption of constant plasma 
properties . 
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Double Ionization 

The presence of doubly charged ions in the discharge chaiiib(;r is 
undesirable for several reasons. Perhaps the most important of these 
is sputtering damage to the interior of the discharge chamber, partic- 
ularly the baffle and screen grid. Therefore, the effect of the 
proposed discharge chamber modification on the double ion concentration 
should be assessed. Theoretical sutdies of the double ionization pro- 
cesses occurring in mercury electron-bombardment thrusters [22] have 
shown that the majority of the double ion production is due to electron 
bombardment of singly charged ions. This permits an approximate 
relationship for the double-to-single ion density ratio to be derived 
based on the assumption that double ions are produced solely from 
single ions. Equating the double ion production and loss rates results 
in the following expression for the double-to-single ion density ratio 


++ 


V/A 


/I V, 


(19) 


The collision frequency v^. is given by Equation (8) with the cross 
secLion o replaced by the cross section for the production of double 
ions from single ions, Comparing Equations (14) and (19) reveals 

the dilemma the ion chamber designer is faced with. In order to achieve 
high propellant utilization efficiencies the discharge chamber charac- 
teristic length V/A should be large, but according to Equation (19) this 
will result in a large double ion concentration. 
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Fortunately, the modification proposed in Figure 7b results in a 
predicted increase in the propellant utilization and a decrease in both 
the discharge power loss and double-to-single ion density ratio for the 
proper choices of the cylinder radius and length. This occurs since 
the open area to neutral propellant loss remains constant (the 
physical open area of the grids remains fixed), while the effective 
open area for ion extraction increases (ion production region 
expanded radially). In other words, the characteristic length V/A can 
be reduced, causing a reduction in the double-to-single ion density 
ratio. At the same time the quantity A^V/A can be increased, producing 
an increase in the propellant utilization; while A/A^ can be reduced, 
producing a decrease in the discharge power loss. As an example, one 
might assume an initial primary electron region length and radius of 
9 cm and 4 cm. If the primary electron region length were reduced to 
4 cm and the radius increased to 6 cm, one could expect a 7 percent 
increase in propellant utilization, and a reduction in discharge 
power and double-to-single ion density ratio of 48 and 14 percent, 
respectively. 


III. APPARATUS 


The discharge chamber used in this invest icp is iiiiisirai i i 
Figure 9 along with the SERT II thruster disc.li.n < ii.iinlier. The mai.M 
difference between the two configurations is the additional soft-iron 
pole piece located between the cathode and anode pole pien^s. This 
modification produces the cusped magnetic field (CMI) ijt.ometr.y 11 Im 
trated in the figure and results in an increase in the neutral resideru 
length throughout most of the chamber. The SERT II thruster permaneni 
magnets were replaced with electromagnets, and current to the magnet 
windings is provided by separate power supplies enabling independent 
control of the fields produced in the regions between the pole pieces. 
The original soft- iron thruster end plate was replaced with an aluminum 
plate, and the upstream electromagnets were designed to provide a 
permeable path between the cathode and center pole pieces. These 
features insure that the minimum path length through the plasma in the 
upstream end of the chamber is between the pole pieces and results in 
a predominantly radial magnetic field in this region. 

The cusped magnetic field geometry is produced by maintaining the 
cathode and anode pole pieces at the same magnetic polarity. Reversing 
the current direction through either set of electromagnets changes the 
polarity of this set and results in a divergent magnetic field geometry. 
This unique feature allows radical changes in the magnetic field 
geometry to be made during thruster operation by simply reversing a 
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switch. The purpose of the divergent magnetic field geometry was not 
to duplicate the SERT II magnetic field, but to allow the eff^-cts of 
two radically different magnetic field geometries to be evaluated with 
the same discharge chamber and identical operating conditions. An iron 
filings map illustrating the field geometries obtained with this 

arrangement is presented in Figure 10. 

The length of the cylindrical anode located in the downstream end 
of the chamber was reduced to permit the installation of the center 
magnet pole piece, and a disc-shaped anode was installed in the upstream 
end of the chamber. This electrode is constructed of stainless steel 
screen which provides a uniform distribution of the neutral propellant 
flow to the discharge chamber. The axial position of the upstream anode 
can be varied during thruster operation, and this capability permits 
some operator control of the plasma properties and thruster performance 
by altering the plasma confinement at the upstream boundary of the 
primary electron region. 

The separation distance L between the center and anode pole pieces 
determines the length-to-diameter ratio of the primary electron region, 
and this parameter is easily varied by replacing the downstream section 
of the chamber. The ability to vary the primary electron region length, 
magnetic field geometry, and anode position makes the discharge chamber 
designed for this investigation a unique and flexible research tool. 

Some additional features include the following: (1) a variable 

magnetic baffle [23] which is useful in regulating the conductivity of 
the plasma in the annular region between the cathode pole piece and 
baffle, (2) individual vaporizers which permit independent control of 
the main and hollow cathode flow rates, and (3) large orifice diameter 



(0.5 mm) main and neutralizer cathodes which facilitate thruster 
operation at high beam current conditions. 

The accelerator system is a high perveance dished type which pro- 
vides good thermal stability and high beam current capacity at the 
specific impulse level desired for North-South stationkeeping of 
geosynchronous satellites. The screen electrode has a 76 percent open- 
area fraction to minimize ion recombination losses, and the accelerator 
electrode has a 54 percent open-area fraction to reduce neutral propel- 
lant loss. The accelerator grid is compeni.ated to reduce the beam 
divergence and resultant thrust loss caused by the spherical shape of 
the accelerator system. 

Mass flow rates are determined by timing the drop of the liquid 
mercury column in a precision bore glass tube, and these measurements 
are considered accurate to within ±3 percent. Currents and voltages 
used in performance calculations are measured using digital voltmeters, 
and all other electrical measurements are made with moving-coil meters. 
The accuracy of the digital meters is il percent of full scale, while 
that of the panel meters is +2 percent of full scale. 

Plasma diagnostic information is obtained by the use of a movable 
Langmuir probe [24]. The collector is 0.074 cm diameter tantalum wire 
which extends 0.127 cm beyond an aluminum oxide insulator. A sliding 
protective sleeve is used to shield the probe and insulator from 
sputtered material when not in use. 

A Faraday cup probe [24] is used to measure the ion beam current 
density profiles at a distance of 6 mm downstream of the accele’'ator 
()rid. The double and single ion components of the curren density pro- 
file are obtained by use of a collimating f: x [5 momentum analyzer [25]. 
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A precision stepper motor is used to pitch the analyzer to within 
±0.05° of the desired angle, and an electrically driven worm drive 
assembly is employed to sweep the analyzer through the ion beam. The 
collector current is measured by a Model 416 Keithley picoammeter 
having an accuracy of ±3 percent of full scale. 

The thruster was operated in the Ion Propulsion Research Facility 
at Colorado State University. The vacuum tank is a stainless steel 
cylinder 1.2 m in diameter and 4.6 m long and has a cryogenic target 
and liner. High-vacuum pumping is accomplished with an oil diffusion 
pump, and the tank liner, target, and diffusion pump baffles are cooled 
by circulating liquid nitrogen. Curing thruster operation the tank 
pressure is in the 10" torr range. 


IV. PROCEDURE 


Thruster operating conditions were selected on the basis of using 
a 15 cm diameter mercury thruster to satisfy the North-South station- 
keeping requirements of a 700 kg geosynchronous satellite [26]. The 
specific impulse and thrust requirements of this mission dictate a 
600 mA beam current with screen and accelerator grid potentials of 
+1000 V and -500 V, respectively [27]. Typical thruster operation at 
the performance curve knee results in a propellant utilization effi- 
ciency of about 80-85 percent. This range of utilization and the 
desired 600 mA beam current results in a nominal propellant flow rate 
of 730 mA which was the standard flow used in this investigation. 

A 37 V anode potential was chosen as the nominal value and repre- 
sents a compromise between ionization efficiency and thruster lifetime 
considerations. The anode power supply was operated in the current- 
limited mode, and the desired anode potential was achieved by varying 
the plasma conductivity in the baffle aperture region. The variation 
could be achieved either by altering the division between the cathode 
and main flow rates or by use of the variable magnetic baffle. In 
general it was found that a slight performance gain could be realized 
by using the cathode flow rate to achieve the desired voltage level, 
and the magnetic baffle was used primarily to stabilize thruster con- 
ditions immediately after startup. The anode potential could then be 
maintained by simultaneously reducing both the baffle magnet current 
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and cathode flow rate. The low cathode flow condition which was 
reached when the baffle circuit was de-energized is apparently respon- 
sible for the enhanced performance and has been observed by other 
investigators [20]. 

Four discharge chamber lengths were evaluated during this investi- 
gation and correspond to 100, 75, 63, and 38 percent of the SERT II 
thruster length. For all but the short discharge chamber, data were 
obtained for both the cusped and divergent magnetic field geometries 
with the upstream anode located in two different axial positions. 

After the initial starting transient had decayed, the propellant 
flow rates were monitored and vaporizer currents varied until the 
desired flow rate was reached and three consecutive total flow readings 
were within ±1 percent of each other. As soon as the desired steady- 
state conditions were reached, data were recorded in the following 
sequence; (1) The Faraday cup probe was swept through the ion beam 
and the probe collection current recorded at 5° intervals. These data 
were later punched on cards and used as input to a computer program 
which performed the integration for the beam current and beam flatness 
parameter. The integrated beam current determined in this manner was 
usually in agreement with the measured current to within *2 percent. 

(2) The Langmuir probe was positioned at various radial and axial 
locations within the discharge chamber, and a continuous variation of 
collector current versus voltage was obtained by use of an x-y plotter. 
The soacing used in the probe surveys was 1.3 cm axially and 2.2 cm 
radially. This gave a total of 28, 20, 16, and 4 probe locations for 
the four chamber lengths investigated, The plotted Langmuir probe 
output was converted to digital form and punched on cards by use of an 


Aulotrol digitizer. These data were then used as .iii»ut. to a computer 
routine which performed calculations to yield the plasma properties. 

A numerical procedure for analyzing the Langmuir probe data was 
developed as a part of this investigation and is described by Beattie 
[28]. (3) The ^ x momentum analyzer was swept through the ion beam 

at pitch angles of 0, 5, 10, and 15°. The collector current sensed by 
a pi coammeter, as well as the output of a position- indicating poten- 
tiometer, was recorded by use of an x-y plotter. The current versus 
position curves were subsequently digitized and used as input to a 
computer routine which performs the double integration required to 
obtain the single and double ion beam current components. The integral 
equations used in the program are derived by Wilbur [22]. (4) Dis- 

charge power versus propellant utilization curves were obtained by 
varying the output current of the anode power supply. During the per- 
formance data acquisition the anode potential was allowed to vary in a 
manner determined by the conductivity of the discharge chamber plasma. 
Although this procedure permits rapid data acquisition and gives some 
idea of the performance trends, it should be pointed out that data 
obtained in this manner do not necessarily indicate the optimum 
thruster performance since they were obtained at a nonoptimum arc 
voltage level. 

The use of a sliding protective sleeve on the Langmuir probe 
insulator was found to extend the useful lifetime of the probe by 
almost an order of magnitude. Unshielded probe insulators often became 
coated with sputtered conducting material resulting in an increase in 
probe collection area and erroneous plasma properties. The sleeve 
assembly effectively shielded the probe insulator from sputtered 
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mdterial when the probe was not in use. By rotating the probe collector 
vertically upward the sleeve would slide downward and expose the 
collector to the plasma. After the probe survey was completed the 
probe was rotated downward and the sleeve would slide down the probe 
insulator and enclose the collector. The probe assembly was stored in 
a low-density region of the discharge chamber which further reduced the 
possibility of contamination and minimized any disturbing effects the 
probe might have on the plasma. 


V. EXPERIMENTAL RESULTS 


The experimental results obtained with the cusped and divergent 
magnetic field geometries are presented in the first and second sections 
of this chapter. The third section presents the results for a radial 
magnetic field geometry which represents the limiting case of zero 
length-to-diameter ratio for the discharge chamber of both the cusped 
and divergent field geometries. The final section presents a compari- 
son between the optimum cusped magnetic field discharge chamber and 
the SERT II thruster. 

Cusped Magnetic Field Thruster 

Beam current density profiles are presented in Figure 11 for the 
three chamber length-to-diameter ratios investigated. Two profiles are 
presented for each chamber length; one corresponds to the upstream 
position* of the anode which resulted in optimum thruster performance, 
and the other corresponds to a downstream, or nonoptimum, anode position 
which degraded thruster performance. The beam flatness parameter was 
inrreased slightly with the anode in the nonoptimum position for all 
three discharge chamber lengths. However, changes in the flatness 
parameter were less than 3 percent suggesting the anode position has 
little effect on the beam flatness parameter. The effect of discharge 

*The positions refer to the upstream anode. The downstream anode 
(e.<i.. Figure 9) was connected electrically in parallel with the up- 
stream anode, and its position was fixed. 
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Figure 11c. Cusped magnetic field thruster beam current 
density profiles (L/D = 0.23). 730 niA 
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chamber length reductions was to increase the beam flatness parameter 
some 12 percent as the chamber length was reduced to 63% of the SERT II 
thruster length. Figure 12 presents the beam current density profile 
for the short discharge chamber operated at a throttled propellant flow 
rate of 461 mA. Comparison of this figure with the corresponding pro- 
file of Figure 11c indicates a substantial increase in the beam flatness 
parameter was realized as a result of throttling the propellant flow 
rate. 

Performance curves corresponding to the optimum upstream anode 
position are presented in Figure 13 for the three discharge chamber 
lengths investigated. The discharge power loss was reduced approxi- 
mately 100 eV/ion as a result of decreasing the discharge chamber 
length-to-diameter ratio, and the short and intermediate chambers have 
about the same performance. The propellant utilization has not been 
corrected for the presence of multiply charged ions, and, as a result, 
utilizations of over 100 percent are indicated in Figure 13. The solid 
symbols on the performance curves indicate a 37 V anode potential and 
correspond to the operating conditions at which the Faraday, Langmuir, 
and f X ^ momentum analyzer probe data were recorded. The flagged 
symbols used in Figure 13 correspond to the nonoptimum anode position 
and indicate the severe performance degradation which accompanied the 
downstream movement of the anode. The close proximity of these data 
points indicates the performance is essentially independent of the 
discharge chamber length when the anode is in the downstream position. 

The ratio of the double-to-single ion beam current components is 
presented as a function of discharge chamber length-to-diameter ratio 
in Figure 14. The downstream position of the anode is seen generally 
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Figure 12. Cusped magnetic field thruster beam current density 
profile (throttled flow condition, L/D = 0.28). 
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Figure 13. Cusped magnetic field thruster performance 
comparison. 730 mA propellant flow rate. 
Shaded symbols indicate the performance 
corresponding to 4 A anode current and 37 V 
anode potential. 
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Figure 14. Cusped magnetic field thruster double-to-single 
ion beam current ratio. 730 inA propellant flow 
rate, 4 A anode current, 37 V anode potential. 
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where I /I is the double-to-single ion be<ini current ratio. These 
corrections are necessary since each double ion contribution to the 
measured beam current is 2/2 greater than the single ion contribution 
due to the higher charge and velocity of the double ion. The propel- 
lant utilization and thrust correction factors calculated for a double- 
to-single ion current ratio of 5 percent, which is typical of Figure 14, 
are 0.976 and 0.986, respectively. These relatively small corrections 
correspond to a 37 V anode potential which is indicated by the use of 
shaded symbols in Figure 13. Propellant utilizations greater than the 
values indicated by the shaded symbols were obtained at higher discharge 
voltage levels which results in an increase in the production rate of 
doubly charged ions. As a result one would expect a larger correction 
to the measured propellant utilization at the higher power levels, and 
this is indicated in Figure 13 by the fact that the measured propellant 
utilization exceeds 100 percent. Attempts to obtain propellant utili- 
zation and thrust correction factors at high power conditions were 
unsuccessful due to the thermal expansion and subsequent shorting of 
the accelerator system caused by prolonged thruster operation at high 
discharge power levels. However, as mentioned earlier, the thruster 
would in general not be operated with an anode potential in excess of 
37 V due to component lifetime considerations. 

Divergent Magnetic Field Thruster 

Figure 15 presents beam current density profiles for the divergent 


magnetic field geometry. The optimum, or upstream, anode position 
generally resulted in a more uniform beam profile and the effect of the 
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Figure 15c. Divergent magnetic field thruster beam 
current density profiles ^'/D = D.?3). 
730 mA propellant flow ra , 4 A anode 
current, 37 V anode potei tial. 



nonoptimum, or downstream, location was relatively i^mall. The Iteain 
flatness parameter decreased when the chamber lennth-to-diaiiieter ratio 
was reduced from 0.53 to 0.3. Further length reduction produced an 
increase in beam profile flatness but not of the same magnitude as the 
original reduction. As a result the highest beam flatness parameter 
was realized with the longest discharge chamber; just the opposite of 
the results obtained with the cusped magnetic field geometry. With the 
exception of the intermediate discharge chamber length, the effect of 
moving the upstream anode to a nonoptimum position was to reduce the 
beam flatness parameter. However, as with the cusped field geometry, 
the effect of anode position was small enough that the flatness param- 
eter can be considered essentially independent of anode position. 

Performance data for the divergent magnetic field thruster are 
presented in Figure 16. The shaded symools represent the 37 V anode 
potential and correspond to about the same power input to the discharge 
chamber. The effect of reducing the discharge chamber length was to 
increase the ion beam current which results in a reduction in the dis- 
charge power loss and an increase in propellant utilization. This 
effect was also realized with the cusped magnetic field geometry. The 
flagged symbols in Figure 16 indicate the performance level with the 
anode in the downstream position and demonstrate the substantial per- 
formance degradation which the nonoptimum anode location produced. A 
comparison of the effects of anode position on thruster performance for 
both the cusped and divergent field geometries shows the diver()cnt field 
thruster was affected to a lesser degree than the cusped field thruster. 
The performance data of Figures 13 and 16 for the nonoptimum atiode 
position suggest the chamber length variation had a greater ('ffect on 
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Figure 16. Divergent magnetic field thruster performance 
comparison. 730 mA propellant flow rate. 
Shaded symbols indicate the performance 
corresponding to 4 A anode current and 37 V 
anode potential. 
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the divergent field geometry than on l.he (n'.|)ed ll('ld geoiiKM.ry . Ihe 
performance data also indicate the baseline disc.liarge power loss for <i 
given chamber length is approximately the same for either magnetic 
field geometry. In addition, the performance curves for the divergf>nt 
magnetic field geometry have a more pronounced "knee" than the corres- 
ponding data for the cusped magnetic field geometry. 

The ratios of the double-to-single ion beam currents are presented 
in Figure 17 for the divergent magnetic field geometry. The effect of 
reducing the discharge chamber length is to increase the double-to- 
single ion current ratio, and the effect is more pronounced with the 
anode in the upstream position. The effect of the downstream anode 
position was to substantially lower the double-to-single ion current 
ratio. Comparison of Figures 14 and 17 indicates the cusped magnetic 
field discharge chamber geometry generally has a lower double-to-single 
ion density ratio with the anode in the upstream location. With the 
downstream anode position the divergent magnetic field geometry has a 
slightly lower double-to-single ion density ratio. 

Radial Magnetic Field Thruster 

The limiting case of a zero length-to-diameter ratio discharge 
chamber was investigated for both the cusped and divergent magnetic 
field geometries by removing the downstream section of the discharge 
chamber. The resulting chamber geometry is illustrated in Figure 18. 

In this configuration the magnetic field is predciidnanLly radi.il and 
the discharge chamber length is 38 percent, of the SFRT II thruster 
length. The end plate used in this thruster was identical to the SFRT 
II design, and the dished g ids v/ere uncoinf)ensated and hafi a 68 percent 
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Figure 17. Divergent magnetic field thruster double-to-single 
ion beam current ratio. 730 niA propellant flow 
rate, 4 A anode current, 37 V anode potential. 
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open-area fraction. These minor configuration changes are assumed to 

★ 

have little effect on the thruster beam profile and performance level, 
and the results obtained with this thruster will be compared with the 
results presented in the proceeding sections. 

The beam current density profile of the radial field discharge 
chamber geometry is presented in Figure 19. The profile is highly 
nonuniform (F = 0.23) and resembles the profile of the axial field 
thruster illustrated in Figure 5a. Performance data for the radial 
field thruster geometry are presented in Figure 20 and indicate a 
serious performance degradation resulted from reducing the chamber 
length-to-diameter ratio to zero. Although the baseline discharge loss 
is comparable to the cusped and divergent magnetic field thrusters, the 
propellant utilization at the knee of the performance curve is only 
about 40 percent. 

Comparison With The SERT II Thruster 

Nondimensional ion beam current density profiles are presented in 
Figure 21 for the SERT II and cusped magnetic field thrusters. The 
integrated beam current is the same for each profile, and when normal- 
ized in this manner the reciprocal of the peak normalized current 
density is numerically equal to the beam flatness parameter. The ideal 
current density profile has a normalized current density of unity and 
is shown in Figure 21 for comparison. The normalized profiles indicate 
a dramatic improvement over the SERT U thruster beam profile uniformity 

*The cusped magnetic field ihruster was operated with these same 
minor changes, and the effects on beam profile flatness and performance 
level were quite small. 
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Figure 19. Radial magnetic, field thruster beam current 
density profile. 492 mA propellant flow 
rate, 3.3 A anode current. 37 V anode 
potential . 
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Compdrison of th6 SERT II and cuspod nidgnotic fiold thrustor 
beam current density profiles. 730 mA propellant flow rate, 
4 A anode current, 37 V anode potential. 


was accomplished with the cusped magnetic field dischanio chamber 
design. The beam flatness parameter of the cusped field geoim-iry is 
about 30 percent greater than the SERT II thruster value for t.lie same 
discharge chamber length. The flatness parameter of the shorter 
chamber geometries is about 40 percent higher than the SERT II thrust, 
value. 

A performance comparison between the SERT II and the cusped field 
thrusters is presented in Figure 22. The cusped field data correspond 
to a discharge chamber length-to-diameter ratio of 0.3 and indicate a 
slightly higher baseline discharge loss than the SERT II thruster. The 
performance at the "knee" of the curves is nearly identical for both 
thrusters. 

A summary of the thruster operating conditions, performance level, 
beam flatness, and double-to-single ion current ratio is presented in 
Table 1 for the SERT II and optimum cusped field thrusters. These data 
indicate the cusped magnetic field thruster offers substantial improve- 
ments in the uniformity of the beam and double ion current density 
profiles. The increase in the flatness parameters of these profiles is 
about 40 percent. In addition, the double-to-single ion current ratio 
in the cusped field geometry is approximately 36 percent less than the 
corresponding SERT II thruster value. The performance degradation 
accompanying these improvements is relatively small; approximately 7 
percent increase in discharge power loss and a 4 percent reduction in 
propellant utilization. 


DISCHARGE LOSSES (eV/ion) 





TABLE 1 


COMPARISON OF THE SERT II THRUSTER AND OPTIMUM 
CUSPED FIELD THRUSTER OPERATING CONDITIONS 
AND PERFORMANCE PARAMETERS 



SERT II 
THRUSTER 

CUSPED FIELD 
THRUSTER (L/D=0.3) 

Operatinq Conditions 
Anode Potential (V) 

38 

4.1 

36.7 

Anode Current (A) 

4 .3 

^ o c: 

Beam Current (mA) 

654 


Propellant Flow Rate (mA) 

725 

730 

Performance Parameters 


1 

Discharge Power Loss (eV/ion) 

239 

256 

Corrected Propellant Utilization (/o) 

86.8 

8.1 

83 «b 

t o 

Double-to-Single Ion Current Ratio (T) 

b .2 

Beam Flatness Parameter 

0.49 

0.68 
r\ on 

Double-Ion Beam Flatness Parameter 

j 0.27 

0.38 


VI. DISCUSSION OF THL EXPERIMENTAL RESULTS 


A discussion of the improvement in the ion beam uniformity which 
resulted from the cusped magnetic field modification is presented in 
the first section of this chapter. The next section describes a pro- 
cedure which was developed for use in defining the ion production 
region of an arbitrary discharge chamber configuration. The following 
section presents the volume-averaged pT^sma properties for the ion 
production region defined in this manner. These properties are then 
used to explain many of the performance characteristics presented in 
Chapter V. A comparison between the experimental results and the 
predictions based on the performance and double ionization considera- 
tions of Chapter II is presented in the next part of the chapter. The 
final section discusses the beam profile shape and performance charac- 
teristics of the radial field thruster configuration. 

Beam Profile 

The cusped magnetic field discharge chamber modification was quite 
successful in improving the uniformity of the SERT II thruster beam 
current density profile. For the same discharge chamber length, the 
cusped magnetic field configuration resulted in an increase in the 
SERT U thruster beam flatness parameter of about 30 percent (from 0.49 
to 0.64). Reducing .he discharge chamber length resulted in further 
improvements in the beam profile uniformity, and with the chamber lenqth 
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rcducod to 63 porcont of the SIRT II tlinistor v.iluc flic (lafncss |iaram- 
eter was about 40 porcont hiijhor. Th(> iniprovoinciit in the ion l.oam 
flatness is due to the oxistenco of a radially uniform ion density 
profile near the screen grid. This is illustrated in Figure 23 which 
shows the ion density profiles for the three different discharge chambei 
lengths. The improvement in the ion beam current density uniformity 
which accompanied the discharge chamber length reductions can be ex- 
plained by considering the magnetic field variations produced by these 
reductions. Figure 24 shows the radial variation in the axial component 
of the magnetic field at a point midway between the center and anode 
pole pieces. These measurements indicate a substantial reduction in the 
magnetic field fringing was realized with the shorter discharge cham- 
bers, and this produces a relatively field-free region over a large 
fraction of the discharge chamber volume. The absence of a significant 
magnetic field eliminates the need for a radial plasma density gradient 
and results in much more uniform ion density and beam current density 

profiles. 


Ion Production Region 

In divergent field thrusters the region of the chamber where the 
bulk of the ion production occurs has been defined by considering an 
iron filings map of the magnetic field configuration. Using this pro- 
cedure the upstream boundary of the ion production region can be 
approximated by the "critical" maonetic field line which is defined as 
the innermost field line which intercepts the cy1 indrii <il anode. In 
the cusped field tlirustor discliarof' cli.imbc'r t ho dilinition o1 the 
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iiiii(|iicl. poll' |ii('c:c (itiiHlf. liri'. I idh )ir(iin|il (•(! ttii' I iip- 

of a (|u<in I, i l.rif i VO prof cfluro loi doliniii'i l.lio ion pi(iil'K i ion volnmo 
iti an arbil.rary (fi oUiarpo (.liaiiihor (lo^il|ll. llio iiioMioil dovi- loprd lor 
(Jcrliniiuj tlio homidary ol l.lio ion prodtKl.ion volmiio (,oii‘. i '.I.'. ol lln' 
follovniv) stops; (1) Tho iflasiiia proportios, ,is flotonniiiod by l.inMninii 
probe measurements, .ire used to calt.ulale l.lii,’ ionization collision lie 
quency at each probe location by use of fquation (H). The c.iluii.ited 
collision frequencies are normalized with respect to the maximum value, 
and contours of the normalized collision frequencies are determined. 

(2) The single and double ion production rates are calculated for the 
entire discharge chamber volume. Tho collision frequency contour whicn 
encloses a large fraction (e.g., 95 percent) of the total single or 
double ion production rate is then used as the ion production re()ion 
boundary. Once the boundary has been identified, volume-averaged 
values for the disctiarge chamber plasma properties can be calc.jlated. 
These volume-averaged properties are useful for comparing the various 
discharge chamber configurations and explaining the thruster performanu' 
characteristics. As a part of this investigation, a computer routine 
was developed to perform the analysis described above. The program w.is 
written to be compatible witli the ion production and vol umo-averaging 
programs of Peters and Wilbur [29] and calls these proep'ams a<. sub- 
routines. Input to 1 fie progr.mi consists of the disih.irge ch.imber 
rieometry, Langmuir iirobe positions, plasma proiunM ie:. . and the spciitiod 
fract ion of tlio tot.il ion profhu I ion Lh.it tlu' final .ontonr line <diould 
enclose. The progr.mi output ((Uisi'.ts ol a gr.ipliiial display o1 the 
dis(li.i*'go ch.imber .|eiiia(d ry and (..riisieii IrciM.-tKy .onlians. ,is v/ell as 
t.tie ion |iro(lu( I ion region biaindary. An (’•siniple ol liie .olii'. ion 
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Iho pr.:.. hr-' I 'vs'ifp'l ahov.> f-vr id sitifyinci l-.ho ion production 
i'is:)ion li'iiiiid.i' il'"'’ 'I'Oii ipiilo ir-'inl in lii frilf'‘(iiK' .'ospectsi (1) 
fiiis iiioLl'iod of idt nti fyiiK] th'"' iiiiii'iri..uil. i I'jioii for ion production 
resulioad in iniicii liotior agn'oiiient bofweon calculdtr*d and iiie.'isured 
'|iian f i L i I's , nirh as propoltant utilisation and doublo'to-singlo ion 
density ratio, than had been realized when Llie ion production region 
boundary was identified using iron filings maps, (i^) Ihe collision 
I n ■' |Uoii(.' y (oiituui's iil'antify f undaiiidifnl dilferences liotweon the cuspcd 
and d i vt’V'jent iiktyiieiic fiehi 'jeoiiifd ries . (or exaiii|ile, tfie tontout 
lines lor Llie divergent magnetic field geometry are seen to diverge, 
imu il like flic iiiagnotic field linos. The conCour lines for the cusped 
iiuiyiici ic iicld geuiiieLry, on the oflior hand, cio.irly iiulicato a phenom- 
enon tliaf i*. not iinmediaiol y obvious from considering magnetic field 
1 i no''t Ihe iniiiloiir 1 inc) ari' drawn out towards the center magnet pole 
pieie iiidiiatin'i .i hi'jli energy deir.ily region extendin') radially outward 
I roi" I ho (ciili i "I Mil' d i '.( lilt rge ( iMiiiher. (3) Itie ion production region 
hoiiMil.ii M , sb'.'i;/ iiiiosii Mil' nil 1 o.i ill the deiiMi of the ion 
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Figure 2S 


. Collision frequency contours and ion production reqion 
boundary (divergent magnetic field geometry, l/D =0.3 
Collision frequency contour increment = 0,1, 
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Fiqure 26 


ION PRODUCTION REGION BOUNDARY 


Collision frequency contours and ion production region 
boundary (cusped inaqnetic field qeometry, L/D = 0.30). 
Collision fr(;quency contour increment = 0.1. 
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production region near the chamber boundary which was realized with the 
cusped magnetic field geometry. 

The shape of the contour lines in the cusped magnetic field con- 
figuration result in a T-shaped ion production region. This phenomenon 
is due to the existence of a "magnetic bottle" in the region between 
the cathode and center magnet pole pieces which can be visualized by 
considering the iron filings map of Figure 10a. Electrons injected 
through the baffle aperture tend to follow field lines which terminate 
on the center magnet pole piece. However, as the electrons approach 
the pole many are reflected by the strong magnetic field which exists 
there. As the reflected electrons approach the cathode pole piece many 
are reflected by the strong magnetic field 'n this region. This trap- 
ping effect is particularly evident in Figure 27 which shows the con- 
tours for the cusped field discharge chamber having a length equal to 
the SERT II thruster length. Contours for the same discharge chamber 
but a divergent magnetic field geometry are presented in Figure 28. 
Comparison of Figures 27 and 28 illustrates the dramatic change in the 
ion production region geometry caused by the change in magnetic field 
geometry. 

Although the cusped magnetic field geometry generally resulted in 
a more uniform beam profile than the divergent field geometry, the dis- 
charge chamber having a length-to-diameter ratio of 0.53 was an excep- 
tion to this rule. Figures 27 and 28 indicate why the divergent field 
geometry resulted in a more uniform beam profile than the cusped 
magnetic field geometry for this particular chamber length. The 
divergent field geometry has the region ol high ion production located 






near the ricreen grid, and the ion production region is spread out over 
a greater fraction of the chamber cross section. 

The improvement in the beam profile uniformity which resulted from 
discharge chamber length reductions was discussed earlier in the chapter. 
The reason for the flatter profile was identified as an improvement in 
the uniformity of the plasma density profile near the screen grid, and 
this was attributed to a reduction in the magnetic field penetration 
into the discharge chamber. The collision freguency contours of Figures 
26 and 27 illustrate this effect quite vividly. The contours in the 
downstream end of the cusped field discharge chamber have about the 
same shape as the magnetic field lines in this region, and a comparison 
of Figures 26 and 27 illustrates the radial expansion of the ion produc- 
tion region which accompanied the reduction in magnetic field fringing. 


Plasma Properties 

The thruster operating conditions and volume-averaged plasma 
properties are presented in Table 2 for the twelve basic discharge 
chamber configurations (two magnetic field geometries and two upstream 
anode positions for each of the three discharge chamber lengths investi- 
gated). The thruster input power is seen to be about the same for each 
configuration. The operating condition which varied is the beam 
current, or propellant utilization, which was affected most bv the 
discharge chamber length and upstream anode position. Table 2 also 
presents the single and double ion collision frequencies and the 
normalized densities and production rate fractions for each atomic and 
ionic species. These quantities were calculated using the computer 
routines described by Peters and Wilbur [29]. The normalized densities 


TABLE 2 

THRUSTER OPERATING CONDITIONS, VOLUME-AVERAGED PLASMA PROPERTIES, AND CALCULATED ION PRODUCTION FRACTIONS 
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indicate that about 50 percent of the heavy particles are ground state 
atoms, 15 percent nietastable atoms, and 30 percent resonance atoms. 

The remaining heavy particles are single, metastable, and doubly 
charged ions. The single ion production fractions indicate that abou^ 

30 percent of the single ions are produced from ground state atoms, 30 
percent from metastable atoms, and 40 percent from resonance atoms. 

The double ion production fractions, on the other hand, show that about 
95 percent of the doubly charged ions are produced as a result of 
electron interactions with singly charged ions. The number in paren- 
theses adjacent to the production rate fractions indicates the fraction 
of the production rate due to the primary electron group. 

The plasma properties and collision frequencies presented in 
Table 2 can be used to explain many of the experimental observations 
presented in Chapter V. For example, with the cusped magnetic field 
geometry (upstream anode position) the improvement in performance 
caused by the reduction in discharge chamber length is due to the sub- 
stantial increase in the single ion collision frequency v. The plasma 
properties indicate the electron temperature and energy v/ere increased 
significantly when the discharge chamber length-to-diameter ratio was 
reduced from 0.53 to 0.3, and this more than doubled the value of the 
single ion collision frequency. Further reduction in the discharge 
chamber length produced a small increase in the primary electron energy 
and a corresponding increase in ■. Results similar to those described 
above are also evident when considering the effects of discharge chamber 
length reductions in the divet'gent magnetic field geometry. Here a 
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The higher collision frequency is seen to be due primarily to an in- 
crease in the Maxwellian electron temperature and primary electron 
energy. A further reduction in the discharge chamber length produced 
a slight reduction in electron temperature and energy, and this lowered 
the single ion collision frequency. However, the ion density increased 
due to the reduced loss rate caused by the lower electron temperature, 
and this resulted in a slight increase in the beam current. 

The double ionization collision frequencies listed in Table 2 can 
be used to explain the variation in the double-to-single ion current 
ratio presented in Figures 14 and 17. For example, the chamber length 
reductions are seen to result generally in an increase in the collision 
frequency due to the increase in electron temperature and energy. This 
observation agrees fairly well with the measured double-to-single ion 
current ratio trends presented in Figures 14 and 17. The downstream 
movement of the anode caused a substantial reduction in the collision 
frequency as a result of the reduction in the Maxwellian electron tem- 
perature and primary electron energy. This produced a substantial 
reduction in the double-to-single ion current ratio as illustrated in 
Figures 14 and 17. 

The performance degradation and reduction in the double-to-single 
ion density rati which accompanied the downstream movement of the 
anode was due to a substantial reduction in the collision frequencies 
V and This was caused by a significant drop in the Maxwellian 

electron temperature and a slight reduction in the primary electron 
energy. Figure 29 compares the Maxv/ellian electron temperature and 
primary electron energy for the two different anode positions and 
indicates the significant reduction in the Maxwellian electron 
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A. UPSTREAM ANODE POSITION 



B. DOWNSTREAM ANODE POSITION 


fitjuru 29. tffect of cinode position on Maxwellian electron tenperatiire 
and primary electron energy (cusped magnetic field geoi'\etry. 
L/D - n.?3). 730 mA propellant flow rate, 4 A anode current 

37 V anode potential . 
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temperature which resulted from the movement, of the anode to the down- 
stream position. The primary electron energy v/as also reduced, liut the 
reduction is localized to the region close to the anode. An explanation 
for the reduction in electron temperature which resulted from the down- 
stream anode movement can be given in the following manner. The 
electrons are trapped in a magnetic bottle region located a distance r 
from the upstream anode. Energetic electrons with a cyclotron radius 
greater than r will be collected by the anode and thus the maximum 
electron speed which can exist in this region is given by the cyclotron 
equation 


V 


max 



( 22 ) 


where V_,„ is the maximum electron speed, and qB/m is the electron 
cyclotron frequency. If a "maximum" speed is defined for the Maxwellian 
distribution function such that a large fraction (e.g., 99.99 percent) 
of the electrons have speeds less than this value, a relationship 
between and the electron temperature can be obtained 
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max 



(23) 


where q/m is the electron charge-to-mass ratio, and is the Maxwellian 
electron temperature. Equating ttie expressions for results in 
the following relationship between the electron temperature and the 
separation distance 
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This expression indicates the strong dependence of the electron temper- 
ature on the separation distance r. The crude model leading to this 
result cannot be expected to be valid for all values of r, of course, 
since the expression was derived assuming the anode was the only factor 
which limits the electron temperature, and clearly this is not the case. 
What one would expect is that for large values of r the anode position 
has essentially no effect on electron temperature, but as the anode is 
moved closet to the field lines the temperature is significantly 
affected. This is just what was observed experimentally. The variation 
of beam current with anode position is presented in Figure 30 which 
shows that the anode was moved in the downstream direction for a con- 
siderable distance with essentially no effect on the thruster perform- 
ance. Further downstream movement resulted in a sharp reduction in the 
beam current due to the electron temperature drop. Equation (?4) can 
be used to estimate the anode position where the electron temperature 
and beam current should drop. Figure 29 shows that the electron 
temperature with the anode in the upstream position is about 4.4 eV. 

The magnetic field strength near the inner edge* of the circular up- 
stream anode is about 15 gauss. Substituting <bese values into Eouation 
(24) results in a value of r equal to 1.5 cm. This agrees fairly well 
with the results of Figure 30 which indicates a rapid drop in the beam 
current at about this anode position. With the anode in the downstream 

The itiFK’r edge was used since current measurements using a seg- 
mented iifistream anode indicate about 99 percent of the currerd to this 
electrode is (irawn hv the inner segment. 
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I’orfoniuince 

The analysis of Cfiapter II suijyested the discharye chamber modifi- 
cation illustrated in Figure 7b would result in a more uniform ion beam 
profile with improved performance characteristics and fewer double ions 
as compared with the configuration of Figure 7a. The performance and 
double ion predictions were based on the assumptions of a cylindrical 
ion production region and constant plasma properties, and these pre- 
dictions generally proved to be valid in spite of the fact that the ion 
production region is not cylindrical and the plasma properties varied. 

As a result of the plasma property variations, the plasma ion production 
cost is net a constant. Therefore, the discharge power losses 
cannot be compared with the predictions of Chapter II since the plasma 
ion production cost is unknown. Flowever, the propellant utilization 
and double-to-singl e ion density ratio can be calculated by use of the 
expressions developed in Chapter II and then compared with the measured 
values of these quantities. The c.ilcul ations were made using Iguations 
(14) and (19) which were derived assuming uniform plasma T>rn|)evties 
exist within the ion production reijion. The fact that the plasma pro- 
perties in actual thrusters have a spatial variation was accounted for 
t>y introducing an equivalent unifonm property region which lias the same 
volume as the nonuni form propL'rty region. Average plasma properties 
were cieFine<l for tlie equivaU'iit region sm ii tint i.lie average valors 
yield the saisr ir>n production r ite as olitained l>y i nterirat i rri the inn 
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production rate per unit volume over the nonuni toriii pn)perty region. 

The area of the uniform property recjinn v/as reduced '.uch that ttie inn 
loss rates calculated using the uniform plasma property va’ ’"‘s were 
the same as the loss rates calculated by integrating the measured ion 
flux rate over the boundary of the ion production region. A detailed 
description of the procedure used tc calculate the volume-averaged 
plasma properties and the plasma nonuniformity factors used to reduce 
the area of the ion production region is given by Peters and Wilbur 
[29]. The results of the propellant utilization calculations are pre- 
sented in Figure 31 which shows good agreement between the calculated 
and measured values. The calculated and measured values of the double- 
to-single ion density ratio are presented in Figure 32, and the agree- 
ment is considered reasonable. The scatter in the data is of the same 
order reported by other investigators [29]. 

Radial Field Thruster 

The radial field thruster configuration, which represents the 
limiting case of zero length-to-diameter ratio, was found to have a 
nonuniform beam profile and poor propellant utilization. The reason 
for this is evident by considering Figure 33 which shows the beam 
current density profile and an iron filings map of the magnetic field 
geometry. The iron filings map shows that the field lines leaving the 
baffle aperture region intercept the screen grid at a small fraction of 
the thruster radius. This results in the confinement of the high-energy 
electrons to a relatively small fraction of the chamber cross section 
much I’ke the early axial field thrusters. In fact, the beam profiles 
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Figure 33. Radial magnetic field thruster beam current density 
profile and iron filings map (small baffle). 492 mA 
propellant flow rate, 3.3 A anode current, 37 V 
anode potential. 
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of the two thrusters are quite similar (e.g., Figures !j.i and 33). With 
the high-energy electrons concentre ced near the chamber centerline, the 
bulk of the ionization occurs in this region and neutrals outside this 
central core are not likely to become ionized. This results in a 
highly nonuniform beam profile and poor propellant utilization. Expand- 
ing the ion production region radially by moving the field line-screen 
grid interception point should result in a more uniform beam profile 
and improved propellant utilization. Figure 34 shows the results of 
shifting the field line interception point to about three-fourths of 
the chamber radius. Expanding the primary electron region volume to 
cover a greater fraction of the chamber cross section resulted in a 30 
percent increase in the beam flatness parameter and a 26 percent in- 
crease in maximum propellant utilization. 

Further improvements in the beam profile flatness and maximum pro- 
pellant utilization were realized with the radial field thruster 
configuration by eliminating the field line-screen grid interception. 
This was accomplished by moving the accelerator system downstream as 
illustrated in Figure 35. This modification allowed the high-energy 
electrons access to the entire chamber cross section and resulted in a 
much more uniform beam profile. At the same input power conditions the 
beam flatness parameter and maximum propellant utilization were in- 
creased by a factor of two as a result of expanding the ion production 
region. 
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Figure 35. Radial magnetic field thruster iron filings 
map and the effect of grid location on beam 
current density profile. 750 mA propellant 
flow rate, 5 A anode current, 37 V anode 
potential . 


CONCLUSIONS 


The most significant results of this investigation are as follows: 

1. The cusped magnetic field discharge chamber design resulted in 
a 40 percent increase in the flatness parameters of both the 
ion beam current density profile and the double ion beam pro- 
file of the SERT II thruster. Based on the reductions in the 
maximum single and double ion current density, the improvement 
in beam profile uniformity should result in about a 40 percent 
increase in the accelerator system lifetime of this thruster. 

The anticipated lifetime improvement is due to a reduction in 
the charge-exchange ion erosion of the accelerator grid and 
reduced sputtering of the screen grid caused by the impinge- 
ment of doubly charged ions. 

2. The performance degradation which accompanied the improved 
beam profile uniformity is relatively small; a 7 percent in- 
crease in the discharge power loss and a 4 percent reduction 
in the propellant utilization of the SERT II thruster. 

3. Analysis of discharge chamber performance and double ionization 
phenomena indicated that, along with the cusped magnetic field 
modification to the SERT II thruster, the discharge chamber 
length should be reduced. Chamber length reductions were found 
to result in a discharge pov^er loss reduction of 100 eV/ion, 


and an optimum chamber length was found. Both of these exper- 
imental results are in agreement with the predictions based on 
the analysis. 

4. A Quantitative procedure was developed and used for defining 
the ion production region boundary and was successful in 
verifying the critical field line concept of divergent field 
thrusters. Better agreement between theoretical calculations 
and experimental measurements of propellant utilization and 
double-to-singl e ion density ratio was realized using the 
boundary defined in this manner as opposed to a definition 
based on iron filings maps. 

5. The collision frequency contours indicate a fundamental 
difference between the cusped and divergent magnetic field 
geometries. The cusped magnetic field configuration tends to 
trap energetic electrons in a magnetic bottle region in the 
upstream end of the discharge chamber. This phenomenon gives 
rise to a T-shaped ion production region in contrast to the 
bell shape of the divergent magnetic field configuration. 

6. The axial position of the upstream anode was found to have a 
considerable effect on the thruster performance and double- 
to-single ion density ratio. Plasma property measurements 
indicated the performance degradation and reduction in double- 
to-single ion density ratio associated with a downstream 
movement of the anode was due primarily to a substantial 
reduction in the Maxwellian electron temperature. 

7. The limiting case of a zero length-to-diameter ratio tor both 
the cusped and divergent maancLic field geometries was 
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investigated, and the porfoniianco and beam profile flatness 
were seriously degraded with this geometry. The source of the 
problem was identified as the confinement of energetic elec- 
trons to the region close to the chamber centerline, and this 
was shown to be due to the magnetic field lines intercepting 
the screen grid. These observations identified an important 
design consideration. Discharge chambers should be designed 
so that the magnetic field lines originating in the baffle 
aperture region do not intercept the screen grid. 
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